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An uncertainty principle underlying the pinwheel structure in the 

primary visual cortex 

D.Barbierit G.Citti | G.Sanguinetti f A.Sarti § ^ 



1 Introduction 



The visual information in VI is processed by an array of modules called orientation preference columns [T] . 

In some species including humans, orientation columns are radially arranged around singular points like the 

^ spokes of a wheel, that are called pinwheels. The pinwheel structure has been observed first with optical 

^^ imaging techniques [2] and more recently by in vivo two-photon imaging proving their organization with 

OO single cell precision [3]. Many hypotheses of morphogenesis of pinwheels during visual development have 

^^ been proposed e.g. [4] [5] and many reviewed in [6] [7], but at the present the problem is still open [8 . In this 

r"| research we do not intend to propose another model of pinwheel formation, but instead provide evidence 

Q^ that pinwheels are de facto optimal distributions for coding at the best angular position and momentum. In 

the last years many authors have recognized that the functional architecture of VI is locally invariant with 

respect to the symmetry group of rotations and translations SE{2) ^ ^lOj ^1] [12] [13] [M] llSj . In the present 

^ study we show that the orientation cortical maps measured in [2] to construct pinwheels, can be modeled 

ti as coherent states, i.e. the geometric configurations best localized both in angular position and momentum. 

The theory we adopt is based on the well known uncertainty principle, first proved by Heisenberg in quantum 

^ mechanics [16] , and later extended to many other groups of invariance [17] [18] [19] . This classical principle 

^ states an uncertainty between the two noncommutative quantities of linear position and linear momentum 

^^ and it has been adopted in [20] to model receptive profiles of simple cells as minimizers defined on the retinal 

plane. Here we state a corresponding principle in the cortical geometry with SE{2) symmetry based on the 

uncertainty between the two noncommutative quantities of angular position and angular momentum. By 

computing its minimizers we obtain a model of orientation activity maps in the cortex. As it is well known the 

pinwheels configuration is directly constructed from these activity maps [2 , and we will be able to formally 

^^ reproduce and justify their structure, starting from the group symmetries of the functional architecture of 

i-H the visual cortex. The primary visual cortex is then modeled as an integrated system in which the set of 

k/ simple cells implements the SE(2) group, the horizontal connectivity implements its Lie algebra and the 

. fH pinwheels implement its minimal uncertainty states. 

^ 2 The functional architecture of VI 

2.1 Simple cells as elements of a phase space 

Receptive fields of simple cells of the area VI are classically represented as Gabor-like filters able to ac- 
complish an analysis both in space and frequency [20] [21] [22] . Each filter is identified by four parameters 
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{q.i->Q2-, \p\-,^) where q = (^1,^2) are coordinates in the 2D cortical layer C, \p\ and are the engrafted vari- 
ables in the sense of Hubel [23^. The coordinate \p\ codes the modulus of frequency and 6 the orientation 
preference, so that they can be interpreted as a choice of polar coordinates in the frequency variables. Since 
the variables of position and frequency are constitutive of the phase space, the set of simple cells is classically 
labeled by parameters of a 4D phase space which represents the cotangent space T*(C) of the 2D cortical 
layer. We also explicitly note that Gabor filters are complex valued, so that the cortical signals related to 
this cells are also complex. 

2.2 Horizontal connectivity and integral curves in SE{2) 

Simple cells are connected by means of the horizontal connectivity with its typical anisotropic pattern 
described in detail in a number of experiments (e.g. [25]). It has been shown that the functional architecture 
of cortico-cortical connections among simple cells strongly depends on position and orientation preference, 
while it presents a weaker dependence on the frequency. Indeed, while studying the connectivity, very often in 
literature only the three parameters (^1,^2, ^) are considered, while the modulus of frequency \p\ is discarded 
[9] [12] [To] [15] . Moreover, the functional architecture is invariant for change of polarity of simple cells so that 
the variable is defined on [0,7r) = P^, instead of [0,27r). Hence the polar coordinates will be expressed 
as (pi,P2) = |p|(— sin(2^), cos(2^)). These parameters are the coordinates of the reduced cotangent space 
T^(C) which can be identified with the elements of the group SE{2) of translations and rotations of R^, 
representing the group of symmetries of the functional architecture of VI. 

The group acts on a vector (xi, X2, (/:') applied in the origin by rotating it of an angle and by translating 
its application point of (^1,^2)' 

(xi,X2)i — >L^q^^q^^e){xi,X2)=r2e{^\^{ M and(/:)i — > {^ - 9) (1) 

where 

_ /cos(2l9) - sin(2l9) 

^^^ ~ V sin(2(9) cos(2(9) 

In particular L(^q^^q^^0^ is the 2D projection of the group action. 

The tangent space of SE{2) has dimension three and it is generated by the infinitesimal transformations 
along two orthogonal translations and the rotation of the group. However, in order to take into account 
the anisotropy of the connectivity, in [12 and [26 just two generators instead of three were considered. In 
particular in [12 it was proposed to choose as generators the differential operators Xi and X2 representing 
only one infinitesimal translation and the infinitesimal rotation. Projecting them on the 2D cortical layer 
through the local action L defined in (IT]) we obtain 



Xi/(xi,X2) = -^— 

dqi 



/(^(gi,0,0)(^l,^2)) =da:J{xi,X2) 



and 



/(^(0,0,^)(^1,^2)) = {xid^^ -X2dxJf{xi,X2). 


These derivations are left invariant with respect to the action L and are simply the directional derivatives in 
the direction respectively Xi = (1,0)-^ and X2 = (^2, —xi)^. The horizontal connectivity is then modeled 
as the set of integral curves 7(5) = {xi{s),X2{s)) of the vector fields Xi + kX2 having starting point at the 
origin of the local coordinates and varying the parameter /c in R (see Fig. |l(b)[ ): 

7'(5) = (Xi+fcX2)(7(5)), 7(0) =0. 

Computing the left invariant derivatives at the point amounts to choose local coordinates around a fixed 
point. This is why the fan is centered in the origin and oriented along the xi-axis in local coordinates. 
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Figure 1: Functional architecture of VI. 



(a) Set of receptive profiles of simple cells with constant 

(^1,^2,^). The 



frequency modulus |p|as elements of the rotation translation group indexed by the triple g 
parameters (^1,^2) code the cell position in the cortical layer and 6 is the engrafted variable of orientation 
preference. The entire set of (^i,g'2,^) constitutes the reduced phase space where the frequency modulus \p\ 
is missed, (b) A fan of integral curves of the vector fields Xi + kX2 by varying the parameter k (left). The 
curves model the connectivity between cells on the cortical plane and are implemented by neural horizontal 
connectivity. The model is in good agreement with the association fields of Field, Hayes and Hess [24 (right), 
(c) The entire set of fan modeling the functional architecture is obtained by applying the group action to 
the fan 1(b)] 



However this choice is not restrictive, since we can always apply the action L^^^ g^, 6*) (^i('5)5^2 (<§)), and 
recover global coordinates, in which the fan has an arbitrary origin q and direction 0. Hence, in the sequel, 
we will always make this choice of local coordinates. 

The cortical activity / = /(xi, X2) defined on the 2D cortical layer is propagated along the horizontal long 
range connectivity, modeled by the 7-curves so that the gradient of the cortical activity will be computed as 



directional derivative of / along the vector fields Xi and X2. 

2.3 The operators in the reduced Fourier domain 

The reduced phase space T^{C) contains all the position variables (xi, X2), but only the angular component 
(p ^ P'^ of the frequency variable. Hence in the Fourier space it supports only functions defined on a semicircle 
of radius Q fixed, and periodical of period tt as depicted in Fig 2(a) The operators acting on these functions 
are the Fourier transforms of Xi and X2 : 



X-^ = —rtsm{2(p)^ X2 = 'id^ 



(2) 



Let us note that the fields Xi and X2 correspond respectively to angular position and angular momentum 
of the TT periodic quantum pendulum [16 . The parameter Vt is specified by the physics of the problem and 
in our case it will express the reciprocal of the correlation length of pinwheels in accordance with [27^ . 




(a) 




Figure 2: (a) k general function u{(p) in the reduced Fourier plane, where only a circle of radius Vt is 



considered parametrized by the angular variable if. The circle is a doubling of the vertical [0, tt] segments 
in fig |l(a)[ u{lp) is defined on this circle with period tt and then its graph is repeated twice. The height of 
the graph represents \u{lp)\ and the color arg('u((/:))). |(b)| The noncommutativity of the group is visualized 
observing that the final point achieved by a translation Xi followed by a rotation X2 is different from the 
one constructed by a rotation followed by a translation. 



3 An uncertainty principle on the functional geometry 

3.1 The noncommutativity of vector fields 

While studying the relation between angular position and momentum, both in the real and Fourier plane, 
we have to take into account that the composition of a translation and a rotation is not commutative: 



r2e 



+ r2Q 



The difference for d sufficiently small is an increment in the direction orthogonal to q (see Fig. |2(b)[ ). This 
property implies the non commutativity of left invariant derivatives X\ and X2. In particular if we compute 
second mixed derivatives we have that X2Xif{xi,X2) 7^ XiX2f{xi^X2). We will call commutator of Xi and 
X2 the difference X^f = X2X1/ — X1X2/, which is a measure of the non commutativity of the space. The 
same property can be transferred in the transformed fields Xi and X2 whose commutator is 



X. 



X2X1 — X1X2 



) = Qcos{2(p). 



3.2 The uncertainty principle 

This noncommutative framework has important analogies with the canonical quantum mechanical case, in 
which the two non commuting operators are position and momentum. In that specific situation the variance 
of the position and momentum operators on functions have a lower bound that has been formalized in 
the well known Heisenberg uncertainty principle. This principle has been formally extended (see [17] and 
references therein) to general non commuting vector fields. Hence we can express here a similar uncertainty 
principle in terms of the vector fields Xi and X2 



AfXi A;X2 > i 



< Xs >; 



(3) 



where ArX is the standard deviation of the operator X and | < X > ? | is the mean value of the operator X 

when evaluated on the function /. For | < X3 > ? | fixed we recognize that, if / is concentrated in the Fourier 
domain so that its variance is small, then the variance of the angular momentum has to be high and viceversa. 
Note that the inequality expresses the uncertainty principle for the tt periodic quantum pendulum between 
the two observable angular position and angular momentum. For cp ^ 0^ the fields Xi and X2 reduce to the 
classical position and momentum operators of the harmonic oscillator and the inequality expresses the well 
known Heisenberg uncertainty principle. 

4 Coherent states as optimal localization in angular position and 
angular momentum 

4.1 Coherent states in the Fourier space 

The minimizers of the uncertainty principle are called coherent states and correspond to the functions with 
the optimal localization in angular position and angular momentum. Let us first compute the reference 
minimal uncertainty "ground" state Uq{(p) by making use of the well-known equation (see [17] ) 

X2U^{^) = 2iXX,u^{^), (4) 

where the scaling parameter A represents the frequency of the pendulum. Its solution can be analytically 
computed and, up to a normalization constant, is the ground state UQ{(f) = e^^^°^(2^). The set of coherent 
states of SE{2) in the Fourier domain depends on the angular parameter and is obtained as a displacement 
of the ground state uq^^p) of an angle 0: 

u^{^)=e^'''^<^^^-^^\ (5) 

Every coherent state represents the most concentrated function for optimal localization in angular position 
and momentum, both taken starting from the reference angle 6. The function u^{(p) is bell-shaped and 



centered around the angle (see Fig. 3(a)). We note that the bigger is A, the sharper is its localization. 



Indeed its standard deviation can be evaluated as A(p ^ }— and its width reduces around (p = diS QX 
becomes big, leading to a delta S{2{(p — 6)) in the limit. 

If u = S{(p) then sm{(p)S{(p) = and the variance of the angular position is null, meaning that it is 
maximally concentrated. Due to the uncertainty inequality doS{ip) has infinite variance i.e. it is maximally 
undetermined. On the other hand, if u is constant, then dou = and the angular momentum is maximally 
concentrated. 

Up to now, the coherent states u^{(p) have been computed using left invariant operators of the group 
SE{2) and hence imposing invariance on the domain of functions. Now we impose invariance also on the 
co-domain of u^{(p) with respect to local phase rotation e*"^'^^ so that each function 'u^((/9)e*"^'^^ is a coherent 
state. These states are solutions of the equation \4\ where the operators Xi and X2 are expressed in terms 



of covariant derivatives instead of the usual derivatives. These coherent states, visualize in Fig 3(c) are 
invariant for any local change of coordinates in the domain and co- domain, modeling a system that does not 
require any knowledge of an external global reference system. 



4.2 Coherent states in the real plane 

The Fourier anti-transform of the function 'u^((/p)e*"^'^^ defined on a circle of radius Q gives the expression 
uq {xi^X2) of the coherent states in the real plane, providing the maps visualized in figure. 
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Figure 3: (a) A coherent state u^{(p) in the reduced Fourier plane for a fixed value of that is a real 
valued function minimizing the uncertainty principle (|3|. |(b)| (top) The standard deviation of the angular 
position A^f^Xi versus the one of angular momentum A^f7X2 for three different values of A: X ^ X = ^ 
(left), X = ^ (center), A <C 2^ (i"ight). For a better visualization the standard deviation is normalized 
AnXi = AfXi/\{Xs)f\^^ i = 1, 2. (bottom) Coherent states in the real plane obtained by anti-transforming 
u^{(p). For A -^ +00 the states u is maximally concentrated in the angular position (orientation) and it results 
u{xi,X2) = e*^^~^^ sin(26')+x2 cos(26')) ^-j^^^^ ^g ^ p^^e harmouic function, whose level sets are rectilinear and glob- 
ally oriented (left). Note that in this case the Xi derivative in global coordinates vanishes i.e. Xiu{xi^X2) = 
meaning that the angular momentum is maximally spread and the angular position is maximally localized. 
In the other case A ^ the state becomes a Bessel function u{xi,X2) = Jq e*^ v ^1+^2 cos(2((^-6'))^^ achieving 
the best concentration in angular momentum, since its level sets are circles with all the possible curvatures 
(right). Note that in this case X2u{xi^X2) = meaning that the angular position is maximally spread and 
the angular momentum is maximally localized. For A between the two limiting conditions the uncertainty 
are distributed between angular position and momentum achieving the equal distribution of uncertainty for 
A = 2^ (center), (c) A covariant coherent state 'U^((/9)e*"^'^^ in the reduced Fourier plane for a fixed value 
of ^, solution of the equation Q with covariant derivatives. The set of states is invariant with respect to 
local change of coordinates both in the domain and co-domain modeling a system that does not require any 
knowledge of an external global reference. The function is visualized with the same coding as in figure [2(a)| 
(d) The Fourier anti-transform of a covariant coherent state (left) as model of a real activity map recorded 
by Bosking et Al [25 (right). 



4.3 From coherent states to pinwheels 

In [2] the pinwheel structure of VI has been reconstructed starting from a set of cortical activity maps 
acquired with optical imaging techniques in response to gratings with different orientations 9. A color image 
has been obtained from gray valued activity maps, associating a color coding representation to preferred 
orientations [28 . Here we apply the same procedure to the set of coherent states (our model of activity 



maps), producing the pinwheel structure visualized at the center of Fig. 4.3 



MBMSi^ 




Figure 4: Gray valued images visualize covariant coherent states on the 2D real plane by varying the angle 
0. The color image in the center has been constructed from gray valued maps, associating a color coding 
representation to preferred orientations, as in ^. 



5 Qualitative and quantitative comparison with orientation corti- 
cal maps 

The image show a qualitative impressive similarity with the original result obtained in [2] [29] [25] among 
many others. Moreover by computing the power spectrum of recorded pinwheels, following [27], a function 
defined on an annulus with radius Vt is obtained, in complete accordance with the result predicted by the 
presented theory, built up on the basis of the geometry of the functional architecture. In our perspective 
the orientation activity maps are modeled by the coherent states, minimizing uncertainty introduced by 
functional geometry. Note that the coherent states are tt periodic as the measured activity maps. To our 
knowledge in the existing models of pinwheels the activity maps are obtained a posteriori, while in our 
approach a geometric model of the activity maps is directly provided and the pinwheels are constructed as 
in the classical experiments. 
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